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SUMMARY 

The Coulomb T matrix is defined as the limiting value of the T matrix for a 
screened Coulomb potential when the screening is slowly turned off. Study of this limit- 
ing process is facilitated by an expansion in spherical harmonics, and the £ = 0 term is 
treated in detail. The amplitude of the T matrix changes rapidly as a function of energy 
in the vicinity of certain critical points. The precise value of the T matrix in these 
regions requires that screening be taken into account, and the appropriate formulas are 
derived. An expression for the T matrix previously derived by Hostler is shown to be 
correct except in the critical regions. 


INTRODUCTION 

A quantity of considerable interest in the formal theory of scattering is the opera- 
tor T, defined symbolically by 


T = V + V V (1) 

E + ie - K - V 

Here K is the Hamiltonian for the system in the absence of interaction, and V is the 
interaction giving rise to the scattering. The total energy of the system is denoted by E; 
the small imaginary term ie , e — 0 + , serves to make the Green*s function 

G = i (2) 

E + ie - K - V 


well defined. 

Generally it is convenient to think of the operator T as a matrix in the momentum 



representation with matrix elements denoted by ( kg | T | k^ ) . It is also useful to introduce 
a complex wave number k, which is related to the total energy by 



where m is the reduced mass; thus the energy dependence of the T matrix may be indi- 
cated explicitly by writing ( kg | T(k) |k^) , or simply T(k). 

For most quantum mechanical systems, an expression for the T matrix in closed 
form cannot be found. The case of a two-particle system with pure Coulomb interaction 
has been studied extensively, however, and recently Hostler (refs. 1 and 2) has been able 
to derive an integral representation for the Coulomb Green's function that can be written 
in terms of hyper geometric functions. From these it is relatively straightforward to ob- 
tain the Coulomb T matrix. 

The resulting expression for T(k), however, has the drawback that it does not ap- 
proach a well-defined limit as k -► ±|kj | or k — ± | kg | and indeed has branch points 
there. This behavior is certainly not correct, for one can show on very general grounds 
that the only singularities of T(k) should be a branch point at k = 0 and simple poles on 
the imaginary k-axis corresponding to the bound state energies of K + V. 

The correct form of the T matrix when k = ±|kjJ is given in reference 3, where a 
similar anomaly in the limiting process | kg | — | k jJ was studied. The difficulty there 
was traced back to the long-range nature of the Coulomb force and disappeared when the 
effects of shielding were included. 

In the present situation the nonphysical branch points are also due to neglect of 
shielding effects. In all cases of physical interest, the scattering of charged particles is 
caused by an interaction which is screened at very large distances. The T matrix may 
therefore properly be regarded as depending on two parameters, e and the screening 
radius R. When the value of T(k) for real k is desired, the correct order of limiting 
processes is e — 0 followed by R — °°. Usually, the ordering is unimportant, but the 
branch points mentioned occur because in Hostler's expression the limit R — °° has been 
(implicitly) taken first. 

The purpose of the present work is to clarify the situation, first, by showing that the 
correct form of T(k) is obtained as k — ±|k^| and the branch points disappear if the ef- 
fects of shielding are taken into account; and, second, by showing that the order of lim- 
iting processes is unimportant except in the vicinity of k = ± | k ^ | and k = ± j kg j or when 
| k^ | = | kg | . (The last case requires rather special treatment and will not be considered 
here; in what follows it is assumed that |kj | ^ j kg [ . ) 
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UNSCREENED COULOMB T MATRIX 


In this section an expression for the unscreened Coulomb T matrix which is essen- 
tially due to Hostler will be presented. His result for the Green f s function for a pure 
Coulomb potential V(r) = Vp/r is 


.2 _ 


2m 


( k 2 G I k l ) 


6(k 2 - k 1 ) ^ mV Q l + M 

k 2 - k 2 ,V (k 2 - k 2 )(k 2 - k 2 )(k- - k- ) 2 


(4) 


where 



(5) 


)2 = , + V -$<?-*) 


k 2 (k 2 - k x ) 2 


V = 


mVQ 

K 2 k 


( 6 ) 


(7) 


From equation (4), the T matrix may be obtained as follows. The operator relation 

G = 1 + T 

E + ie-K E + ie-K E + ie-K 

becomes, in the momentum representation, 


( 8 ) 


( k 2 I G | k i> 


2m 5 (*2 - k l> / 2m\ 2 <W> 

K 2 k 2 - k 2 \K 2 / (k 2 - k 2 ) (k 2 - k 2 ) 


(9) 


If this equation is solved algebraically for the T matrix and equation (4) is used, 


<Wl>— 1“=^ 

2ir (k 2 - k x r 


( 10 ) 
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The integral M may be evaluated in terms of the hyper geo metric function by 

changing to 

s + 1 

as the variable of integration, with the result 


( k 2 T k l> = 


2tt 2 (k 2 - k x ) 2 



2 F i (h W, 1 + W, £ - ± ^- 
\ P ~\ 


- 2 f i ^ 1 + w. 


(ii) 


Considered as a function of k, the right side of equation (11) has simple poles at 

2 2 

irj = -n (n = 1, 2, 3 . . . ) and branch points at p = 1 and p = These latter points 
correspond to k = ± | kj | , k = ± | kg | , k = 0, and k = The detailed behavior of T(k) as 
k — ± | k j | may be determined by analytic continuation of the hypergeometric series and is 
given by 


<k 2 |T(k)| kj) 



( 12 ) 


where 


-7 r < arg^kg - k 2 ) < 0 
7T < arg(k 2 - k 2 ) < 2 tt 

(The corresponding result for k — ±|k 2 | is readily obtained from the symmetry of the T 
matrix with respect to k^ and kg. ) This result is to be contrasted with the exact result 
of reference 3: 


( k2 I T(ki) | kf ) 

2tt 2 


fe - sfT 


( 13 ) 
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Equations (12) and (13) agree in their dependence on k 2 and (k 2 - k^) , but they differ in 
magnitude and phase (Cq and 8 Q are defined in the next section). In fact, the quantity 
in braces does not approach a well-defined limit as k — ±|kj | , and its magnitude may 
vary by as much as a factor of e vr ^ . For a precise definition of the T matrix in this 
energy region, the effects of screening must be considered. 


SCREENED COULOMB T MATRIX 

Ideally, one would like to find the T matrix for some screened Coulomb potential 
and study its behavior as R — °°. Then, as discussed in the INTRODUCTION, the Coul- 
omb T matrix could meaningfully be defined as 

T(k) = T (R) (k) 

or by a suitable limiting procedure when this limit does not exist. In principle, such a 
program could be carried out by means of an expansion in spherical harmonics, although 
it is highly doubtful that the series for T^ could be summed in closed form. One may 
conjecture, however, that the series for T^°°^ is summable and in fact is given by 
Hostler's expression except when k — ±|kj| or k — ±|k 2 |. Some support for this belief 
may be found in the fact that equation (12) has the right angular dependence with only its 
magnitude and phase in doubt. 

To establish the conjecture rigorously requires examining the general term in the 
expansion of T^ and comparing it to the general term in the expansion of Hostler's re- 
sult. Here a less ambitious program will be followed, that of comparing only the It = 0 
terms. A qualitative investigation of the terms for SL > 0, however, will also be made. 


Expansion In Spherical Harmonics 

In order to parallel the treatment of reference 3, delta-function normalization for the 
wave functions will be adopted, and the expansion of the T matrix will be written in the 
form 


2 00 

< k 2 |T| k i> « - Yj < 2 « + « P A • k l> < k 2l T «l k l> < 14 > 

4ir 2 m f=0 


ik 
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where 


< k 2 l T J k l> ="— < / 1 (k 2 | T |ki>P £ ( M )d M ( M =k 2 • k t ) (15) 

K 2 - 1 

and P £ (ju) is a Legendre polynomial of order f . 

The coefficients T^(k) for an arbitrary potential V(r) may be obtained from equa- 
tions (1) and (15) if the corresponding Green*s function is known. This is accomplished 
by making an expansion of the coordinate representation of the Green*s function: 


( r G r') = 


7Z(^) Pj<; - ; ' )<r|G£|r ' : 

n £=o 


After the angular integrations are carried out, 

/ OO 

j £ ( k 2 r ) w ( r )j f ( k i r )r 2 dr 


where W(r) = (2m/li )V(r) and j^(kr) is a spherical Bessel function of order £ . 
To obtain ( r J G^ | r f ) , the operator equation 

(E + ie - K - V)G = 1 

is written in the coordinate representation, from which it follows that 


lj!l r + k 2.«i±i)_ w(r) 

r dr 2 r 2 


The solution to this equation is well known to be 


(rlGjr.) 

a 

r 


(16) 


/ oo OO 

r 2 dr J r' 2 dr» j £ (k 2 r)W(r) { r | G £ | r* ) W(r*)j £ (k 1 r») (17) 


(18) 


(19) 


<r|G £ |r*> F (r )H £ (r^ 

ikrr* 


( 20 ) 
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where r < is the smaller and r^ the larger of r, r*, and where and are the 
regular and irregular solutions of 



, k 2 _ m + 1 ) 

2 


- W(r)j 


u £ (r) = 0 


having the asymptotic forms 


F,(r) 


^(r) 



( 21 ) 


( 22 ) 


The phase shift 6^ is determined by the fact that F^(r) must vanish at r = 0, and the 
normalization is so chosen that the Wronskian of F^ and is equal to ik. 

By combining equations (17) and (20), one may write 

(k 2 |Tfil fc l> ‘-h'h < 23 > 


where 



j £ (k 2 r)W(r)j jf (k 1 r)r 2 


dr 


(24) 


and where 


3 SL = J £ (12) + J £ (21) 


J £ (12) 



j £ (k 2 r)W(r)F £ (r) 



j £ (k 1 r*)W(r , )H £ (r’)r‘ dr* r dr 


(25) 


and J £ (21) is identical to J £ ( 12) except that kj and k 2 are interchanged. 
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Cutoff Coulomb Wave Functions 

The simplest way to screen the Coulomb field is to cut it off at r = R, so that the 
potential is 


V(r) = V Q /r (r < R) 

= 0 (r > R) 


(26) 


The solutions of equation (21) here must be proportional to pure Coulomb functions for 
r < R and to free-particle functions for r > R. The free-particle functions are so 
chosen that the asymptotic forms of equation (22) are obtained for large r. Thus, 


F (r) 

— = N £ (2kr)^e ikr $(jf + 1 + irj, 2SL + 2; -2ikr) 


1 e^hMkr) + e l6£ h (2) (kr) 
kr 2 L ^ 


H (r) 

-i — = N*(2kr) 1 e lkr ^(f + 1 + it], 21 + 2; -2ikr) 
kr 3L 


kr 



(kr) 


(r<R) 
(r >R) 
(r <R) 
(r >R) 




J 



(27) 


(28) 


where <f> and \I> are confluent hypergeometric functions (ref. 4), and h^ and hj^ are 
spherical Hankel functions of the first and second kind of order k. 

The constants and 6^ may be determined by equating logarithmic derivatives of 
at r = R; approximate values for large R are (ref. 3) 

N f ~ C £ (t?) 6 ! ~ a k ~ * l n ( 2kR ) (29) 


where 


C 1 (t?) = e 


-- irrj-io. 
2 x - 


T(f + 1 + irj) 

T(2k + 2) 


(30) 
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e 


(31) 


2l0 £ _ r(£ + 1 + iT7) 

ru + 1 - iv) 

An exact value for the product will be necessary and can be obtained from the 

fact that the Wronskian of and H £ must be equal to ik. After some algebra, one 
finds that 


F,H. - FJHj = 


*N,Njr(2£ + 2) 


2ir(£ + 1 + i? 7 ) 


(-1)* = ik 


which yields 


— TITj+ifJ £ 


N £ N £ = 2 ( _1 ) £+1 1 2 (-l) £+1 e 2 c 


r(2£ + 2) 


Av) 


(32) 


Finally, from the integral representations of $ and ^ , the following expressions 
(valid for r < R) can be derived: 


W(r)F 0 (r) = - 



W(r)H 0 (r) = - 


k2 Nfr 

2tt 


e 1 Tr? r(l - irj) 


/"(+ 1) 

J (t ^ 8 

V oo 


ds 


(33) 


(34) 


The t-contour is a loop encircling the points t = ±1 in the positive sense; the phases are 
such that arg(t ± 1) = 0 when t crosses the real axis to the right of t = 1. The 
s-contour is a “loop** which begins at s = encircles the point s = 1 in the positive 
sense, and returns to s = the phases are such that initially arg(s ± 1) = 0. 


Evaluation of Tq( k) 

For the cutoff Coulomb potential, the integral Iq becomes 
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(35) 


j - 2rjk f R 


sin k 2 r sin k^r — 


k* + k, 


= ifrr ir^ki + Ci C k i k 2 |R ] ' C1 K + kz)R ] 
1 2 1 2 


where Ci(x) is the cosine integral (ref. 3, p. 267). Since the case k^ = kg is not being 
considered here, 


_ 77k 


In 


k l + k 2 


k l k 2 k l - k 2 


+ & 


ft) 


(36) 


To evaluate Jq, it is convenient to express sin(kgr) and sin(k^r*) in terms of ex- 
ponentials and to introduce the integral representations (33) and (34). Then Jq can be 
written as the sum of four integrals of the type 


■«/“(«)"/ etif 


B dt ds 


(37) 


where B = B(12) + B(21) and 


f R -ir(kt+k 9 ) /' ±t ir*(ks-k.,) 

B(12) = - J e J e 1 dr* dr 


,R irtf 


(38) 


167T 3 i 4ffi 


(39) 


Specifically, 


J 0 = — [A(kj,k 2 ) - A(k 1 ,-k 2 ) - A(-kj,k 2 ) + A(-k 1; -k 2 )l 
k l k 2 


( 40 ) 


The integral for B(12) is easily carried out and finally leads to 
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Figure 1. - Integration contours for confluent hypergeometric functions. 


B = a + 5 + + _ j3 +J e i(a-/3)R 

aa(a - /3) aj3 aft - (3) 


where 


a = ks - k. 


/3 = kt + kg 


(41) 


(42) 


and a and j3 are the corresponding quantities with k^ and kg interchanged. Since B 
has no singularities for finite s and t, the contours may be deformed as desired; de- 
forming them as shown in figure 1 is convenient because then the singularities of the in- 
dividual terms in B do not cross the contours as Im(k) — 0. 

The roles of the four terms in B are as follows. Generally the main contribution is 
from the first term, the others being of order l/R. When k — k^, however, the second 
term cancels part of the first and makes a contribution of its own. A similar remark 
holds for the third term when k — kg. It may be shown (see appendix) that the fourth 
term is always negligible provided that k^ kg; this term will not be considered further. 
Thus in an obvious notation 


and, correspondingly, 


A — A^ + Ag + Ag + 



» tt m 

J 0 = J o + J o + J o + 



(43) 


(44) 


The t-integration for A^ is performed by expanding the contour into a circle of in- 
finite radius, which yields 
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( 45 ) 


(t + r\ 

dt _ 2ni 

/s - v + 1> 

r->] 

\t - 1 > 

' a - 13 k 

\s - v - V 


where v = (k^ + kgVk. Then can be written in the form 

A = 2jriQ lim „ / <+1> ^ 1 -V T| - ll 2 J-LZ 

* ^ *L \s + 1/ [_\s - v - 1/ J (ks - k^)(ks 


- k 9 ) 


ds 


( 46 ) 


Actually the integral is convergent even when L = °°, but the limiting process allows A^ 
to be considered as the difference of two otherwise divergent integrals, which may thus 
be treated separately. In the first of these integrals, let 


x = 


_ (L + 1)(L - v - l)(s - l)(s - v + 1) 
(L - 1)(L - v + l)(s + l)(s - v - 1) 


and in the second let 


y - (L + l)(s - 1) 
(L - l)(s + 1) 


Then, after some manipulation, Aj may be written 



where 


a _ (L - 1)(L - v + 1) 

(L + 1)(L - v - 1) 

> 

b = k^J_ 

L + 1 J 


( 48 ) 
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( 49 ) 


The integrals in equation (47) are all of the same general type, being integral represen- 
tations of a function here denoted by D(x): 

D(x) a (e~ 2m] - l) 1 /* ( ° +) t i7? ■ x - dt 

J j 1 - xt 

(50) 

= — 1 — X oF-iUj 1 + ir i; 2 + i7jj x) 

1 + irj 



Therefore 

A x = i 7 7 k L 1 ^ n oo |a ir} [D(a) - D(auv)] - b i?1 [2D(b) - D(bu) - D(bv)]| (51) 

From the analytic continuation of the hypergeometric function, 

D(x) \ ’/'(I + i7) - ’M 1 ) - ln(l - x) |arg(l - x) | < it (52a) 

D (x) £ -Lp^C^rj) 2 ^ 177 - lj 0 < arg^j < 2n ( 52b ) 

where xf/(x) is the logarithmic derivative of the gamma function. The limit L — °° may 
be taken with the help of equation (52a), and after the combination of equation (40) is per- 
formed 
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Note that the first term in J^, which arises from the limiting value of D(a) - 2D(b), is 
just equal to -Iq. 

When k - kp the expression for simplifies because the arguments of D 
approach either zero or infinity. From equations (50) and (52b) it follows that 






'CM* 


2ik, 



(54) 


with the phases of both l/uv and v/u between 0 and 2n. 

The t-integration for Ag is carried out as for and leads to 



(55) 


where 


f(k x ) = e 



e ikRs 
s - kj/k 


ds 


(56) 


The expression for then takes the form 

[%) - f(-k x )] (57) 

IV IVjXVg 

From the way in which the t contour is laid, one can show that both arg(l/v) and 
arg(v) must be approximately between -37t/2 and ir/2. 

The integral in equation (56) is generally of order 1/R, as may be demonstrated by 
integrating by parts. This argument fails, however, when k - k^ because the contour 
must pass between the singularities at s = 1 and s = k^/k, and these singularities 
merge in the limit. The difficulty may be avoided by first enlarging the contour so that 
it also encloses the pole at s = k^/k: 
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The contour r is now such that, for sufficiently small |k - kjJ, the factor (s - k^/k)~^ 
may be expanded about s = 1 and term-by-term integration performed. The resulting 
integrals may be recognized as confluent hypergeometric functions, and to lowest order 


f(k^) = -2iri 



— — *(i?7, 1; -2ikR) 

r(i - iv) 


+ & [(k - 



(59) 


From the derivation it is clear that 0 < arg(l/u) < 2ir, and since k^ is positive, this 
becomes n < arg(l/u) 2 it. For large R, the asymptotic form of may be used to 
write 


f(k^) = -2jri 



3 

— - ITT] 

— (2kR)' 1TJ 

r(l - ir?) 


+ & 


(i) + e '[< k - k i) R ] 


(60) 


Two things are noteworthy here. First, there is the appearance of the logarithmic 
phase factor (2kR) rr> , known to be necessary when k — |kjJ from the work of refer- 
ence 3. Second, the order of magnitude estimates of the neglected terms in equation (60) 
show clearly that this result can be achieved only if the order of limits is k — k^ fol- 
lowed by R — 

Through use of the definitions of Cq(? 7 ) and 6 q, the final result for J^ f when k kj 
may be expressed as 


J o -- 


C Q (r])e 

2ikr 


i6 r 




(61) 


A similar treatment may be made when k — -k-p for all other values of k, Jff = &(1/R). 

Finally, examination of the third term in B (eq. (41)) shows that J*** is identical to J»» 
except that k^ and k .2 are interchanged. Thus the evaluation of Jq to order l/R is 
complete, and Tg(k) may be determined. 

In the general case, where only Jq is nonvanishing for large R, 


( k 2 1 ^0 I k l ) 


_ rjk 

2k l k 2 



(62) 


in the limit R — °°. This equation corresponds exactly to Hostler's result, as will be 
shown in a later section. In particular, there are branch points at k = ±k^ and k = ±kg 
Clearly, however, these branch points do not appear if the limit R — °° is taken last. 
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Consider the case k — k^; here both and J^’ are important, and when equations 
(54) and (61) are combined 


(k 2 Tq kj) 


C Q {ri)e 


i6. 


2ik r 


[ef-l-a 


( 63 ) 


The limiting value of Tq is independent of the way in which k — kj and, in view of the 
restrictions on arg(l/v) and arg(v), may be written 


< k 2 1 T 0 (k l> I k l > 


C 0 (t])e 

2ikr 


i6. 


k 2 " k l 


k 2 + k i 


177 


k 2 + k i 


k 2 ” k l 


iif* 




Vx 


1 


> + & 


for 


( k 2 > k i 

I k 2 < k i 


(64) 


This is identical to the result found in reference 3. 


PROPERTIES OF T f (k) 

Although no detailed investigation of T^ has been made for higher values of £, 
reasonable predictions about what to expect are not difficult to make. For higher £ val- 
ues in the shielded case, j^k^r*) and j^kgr) would be written in terms of spherical 
Hankel functions; this would again lead naturally to four integrals, of which only the 
first would be independent of R. Presumably the other three would be of order l/R ex- 
cept in the vicinity of the branch points, as in the £ = 0 case. It is not unreasonable, 
then, to suppose that the R-independent term is the same one that would be obtained for 
the unscreened potential, and therefore that the screened and unscreened T matrices 
are identical in the limit R — °° save for the exceptional points. 

To strengthen this argument, consider the integral 

rR rR 

ikj^(12) = J co^r) J o^r') dr* dr (65) 

which arises in the calculation of T«. Here 
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^(r) = rj £ (k 2 r)W(r)F £ (r) (r < R) 

« 2 (r) = rj^k^W^H^r) (r < R) 


( 66 ) 


and because r and r* are always less than R, and are pure Coulomb functions. 
Now the integral in equation (65) may be rewritten 

ikJjj;( 12) ~ w ]( r ) /: y^r^dr* dr - Lr^HU” 0) 2 ( r )drj 


/• oo /*o 

X “> <r) / 


coJr^dr* dr (67) 


provided that <x > ^ and u> 2 are given some suitable definition for r > R. For the present 
purpose, it is convenient to require that and u> 2 have the same functional form for 
r > R as for r < R, that is, and u> 2 are proportional to pure Coulomb functions 
multiplied by spherical Bessel functions for all r. 

With this definition, it is clear that the first term in equation (67) is just what one 
would write for the unscreened Coulomb potential. The second term is similar to 
of the preceding section in that it factors into a function of k 2 multiplied by a function of 
kj. Regarding the second factor, one can show that 



cu 2 (r)dr = 




( 68 ) 


This follows from the fact that for large values of r, 


u> 2 (r) 



e ikr 

(2kr) i7J 


(69) 


The rapid oscillations of the exponentials make the integral small unless k — ±k^. In 
like manner, by using the asymptotic forms of c and one may show (see appendix) 
that the third integral in equation (67) is of order l/R unless both k ±k^ and k — ikg, 
which is possible only if | k^ | = | kg | • 

These findings reveal that J^(12) departs from its unscreened value only when con- 
tributions to the integral from large r are important, and these contributions occur only 
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when k — ±k^. Similar conclusions may be drawn for J^(21), the critical condition be- 
coming k — ikg. Thus, 


t(R) 





(70) 


To complete the investigation of T^, the behavior of 1^ as R — «> must be deter- 
mined. The subtraction procedure used above is also effective here, and leads easily to 


4 R> ■ ♦ - 


( k l " k 2) 


R 


From this and equation (70), it follows that 


t ( r ) = + O’ 


( k±k i) 


R 


+ &\ 


( k * k 2) 


R 


+ 0 \ 


( k l ' k 2) 


R 


(71) 


(72) 


which is the desired result. 

It is instructive to examine the way in which T^ changes as the contributions from 
large r come into play. Consider the case k — k^; from equation (69) and the fact that 

~ 2(-l) f+1 e ,r?7//2 e 1CT£ 


and 

8 £ ~ o- £ - n ln(2kR) 

as R — °° } 



co 2 (r)dr = 2kR f 

kl Jr 


-i ( k i -k ) r dr 

e — 

r l+iT? 



This integral may be evaluated from the formula 


(73) 
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*(1, 1 - in; AR) 


(|arg(A) | < 1/2 tt) 


R i? ? f 

*'R 


-Ar dr -AR 
e = e 

l+i?j 


=-T e 

iV . 


= A. e ~ XR $(1, 1 - irj; AR) - (AR) 1?7 T(1 - i7?)| 


(74) 


The lower form of equation (74) is appropriate to the present case, which corresponds to 
R ~ AR — 0. Because of the restriction on arg(A), 
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A \ i?? = 

i(k x - k) 

\2k/ 

2k 


\ n n\* 


< 2tt 


' + <?(k - kj) i t < arg^ 

Then, since <I>(x) -"1 as x — 0, 

J u> 2 (r)dr = -ie^ + ie 27n? C 0 (7]) e l(C7£ + <y|(k _ k^R] + 


(75) 


(76) 


It will now be shown that 


~i 5 0 

/ at^rjdr = -ke 1 ( k 2 1 T^(k) | k) 

" n 


(77) 


that is, the integral is proportional to T^(k) for the special case when k^ is complex 
and equal to k. For the proof, it is necessary to consider a ’‘plane wave” cp _ whose 

k 

wave number is complex: 

-1 CO 


<Mr) = (2ir) 2 Yj i £ ( 2£ + !) p £ (k • r)j £ (kr) 
k f=0 


K 2 k 2 

2m 


= E + ie 


(78) 


19 



When the scattering operator T(k) is applied to this function, the result may be ex- 
pressed as 

T(k)$>_* = V(r)tMr) (79) 

k k 


where 


V/-(r) = (1 + GV)cp~ 
k k 


= 1 +. 


1 

E + ie - K - 



Thus the T matrix may be written 


(80) 


-3/2 / ik.r 

( k 2 |T(k)|k) = (2ir) Je 2 V(r)^(r)dr ( 81 ) 

k 

The reason for introducing cp becomes clear when one looks for a way to deter- 

k 

mine \1/ . Since cp is annihilated by the operator E + ie - K, it follows from equa- 
ls k 

tion (80) that 


(E + ie - K - V)^_ = 0 

k 


(82) 


Thus is an eigenfunction of K + V belonging to the complex energy E + ie, and 
k 

well-known scattering theory techniques lead to the expansion 

OO 

- -3/2 'ST' c - - iS„ F^kr) 

(r) = (2v) / i £ (2f + l)P«(k • r)e £ — (83) 

k /— / *- kr 

1=0 

The angular integration in equation (81) may now be carried out and yields 

k 2 . i6 # f°° F 0 (kr) 2 

(k 2 |T(k)|k) =— — — (2f + l)P £ (k 2 . k)e 1 J j £ (k 2 r)W(r) — " r dr (84) 

4ir 2 m £=0 0 ^ 

When this equation is compared with equation (14) with the cutoff Coulomb potential for 
V(r), the result given in equation (77) is obtained. 
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One final point remains to be made. In equation (77) the upper limit cannot be ex- 
tended to infinity for arbitrary complex k, because the integral will diverge. It can be 
shown, however, that 



(Oj(r)dr 


k l^ 


(k - k x )R - 0 

R — oo 


(85) 


where g^kg, k j) is a finite, well-behaved function provided that kj A kg. Combining all 
these results yields 


(t^ r)d )(C o> 2 (r)dr^ =ik(k 2 |T ( £ R >(k)|k) 


+ ie 27r7? C 0 (7 ] )e l(CTra ° ) g j? (k 2 ,k 1 )^y 7? + <?[(k - k^] + (86) 

i6 jf 

The first term on the right side is R-dependent through the phase factor e , and 
provides the correct definition of as k — kj. The second term is independent of R 
and has a branch point at k = k^; clearly its function is to cancel a similar quantity from 
J^ 00 ). This can be seen quite easily in the £ = 0 case, since J^°°^ = JJj of the preceding 
section. Hence, 


J 0 R) = J 0 - < k 2l T 0 R) l k > - k C ° (li) )\( k 2’ k l> + *[< k - k l> R ] + «-(^) 


(87) 


From equations (85), (77), and (63), it follows that 


go(k 2 , k i) = 


k 1 C 0 (?7) 


2ik r 


l\ iT ? 


0 


Av 


( 88 ) 


and thus, using equation (54) for Jfo gives 

4 R> = -h - < k 2 l T 0 R) I k > + ^[( k " k l> R ] + (89) 

The preceding discussion indicates that J^°°^ may be written as -ij^ plus a term to 
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be cancelled by contributions from large r. Examination of equation (76) reveals that 
this cancellation is equivalent to replacing 


le 2 ">C 


by 


. i6 f 
le 


as a factor in the term to be cancelled. Since 


6 I m 5 0 * (<r £ - °0> 


i(“) i 


the !. dependence of ' is not affected by this replacement, and 


i5g - 2nrj 


T< R > = ' u 1, 'T < “ ) + M 

1 CM) 1 \R 


as (k - k^)R — 0 


(90) 


The importance of this result is that the angular dependence of will be the same as 

that of T^°°) in the limit k — kp and thus from equation (12), 


<k 2 ]T (R >|ki> =-4 C 0 (»)e 


V 0 15 0 ( k 2 - t 2 ) 1 ’ 1 

h- \ - 0 r /„\ e U \ 2 / 




[(k 2 - kp 2 ] 


+ <y\ 


21 1+i V 


^ as (k - kj)R - 0 (91) 


This result agrees with the findings of reference 3 and in the limit yields equation (13). 


CLOSED FORM FOR T (oo) (k) 

In the foregoing section, it was more or less tacitly assumed that T^°°\k) is given 
by Hostler's formula provided that k is off the real axis. A study of Hostler’s deriva- 
tion indicates that this is true, since he begins with the same Green’s function integrals 
except that the upper limit is °° instead of R. One can also verify directly that 
( kg | Tq | kj ) is the same for both by evaluating 
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1 + M 
(k 2 - k x ) 2 


dp 


(M = kj_ • k 2 ) 


(92) 


An expression for 1 + M is given in equation (11), but a more convenient one for integra- 
tion purposes is 


1 + M = 


1 o 2 -l 
4(1 + i?7) p 


2 F 1 2 > 2 + irr > " zFl (*’ 2; 2 + i7?; ^ 2 "^) 


(93) 


which may be obtained from equation (11) by making use of the transformation properties 
of the hypergeometric function. Then, since 


and 


dp k l k 2 p 2 - 1 

(k 2 -k l} 2 P 


— — f 2 F i (l» 2; 2 + W, x) dx = 

1 + irj J 1 \x - 1/ 


where D(x) is defined in equation (50), the integration yields 

rjk 


< k 2 1 T 0 1 k l > = 


2k l k 2 


■fcr) - ”K) 


/!=+ 1 
M=-l 


r?k 


2k l k 2 


D(uv> + D (™) ' D (v) ' Ku) 


(94) 


(95) 


(96) 


This equation is identical to equation (62), which gives the result for ( k 2 |TQ|k^) in the 

shielded case when only the R-independent terms are included, that is, the result for 

( k 2 1 Tq°^ | k i ) . Although it is beyond the scope of this report to examine in detail the 

exact form of t 1°°^ for H >0, the evidence presented here leaves little reason to doubt 
SL 

the equivalence of T' ’ and Hostler’s formula. 
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CONCLUDING REMARKS 


The Coulomb T matrix T(k) is obtained from the T matrix T'^(k) for a screened 
Coulomb potential by the prescription 

T^(k) - T(k) + (R - oo) 

Generally, it makes no difference whether k approaches the real axis before or after 
the limit R — °° is taken, and in this general case 



V 0 1 + M 
2v t 2 (k 2 - kj) 2 


where M is given by Hostler’s mtegral (5). 

When k = ±|kj | or k — ±|kg | , however, the magnitude and phase of this equation 
are incorrect. Here it is important that k be allowed to approach the real axis before 
the limit R — The result for k — ±|k J is 


___ ___ i6 0 -27T77 __ ___ 

<k 2 |T (R) |k ) - e - — u^kglT^I^) 


2 ^ 


V n i6 n (k« - k*) n 

- C^(r/)e 0 + ^[(kikjR] + ffll 

2 ” 2 1r 


The symmetry of ( kg | T | kj } may be invoked to obtain the result for k — ± | kg | . 

The Coulomb T matrix is seen to be independent of screening effects except when 
k — ± | k j | or k — ± | k g | ; at these points it depends directly on R only through the loga- 


i8 r 


rithmic phase factor e . The magnitude of the T matrix becomes discontinuous in 
the limit of zero screening: for example, if k is on the real axis in the vicinity of kp 


<k 2 T k l> = 


2tt 2 (kg - k^ 2 
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where 


\ 




CqC 7 ?) 

^ > k 

1 1 

k^ = k 


kjCk 

C 0 to) 

[k 2 > k. 


i k 2 < k. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, February 5, 1965. 



APPENDIX - ORDER OF MAGNITUDE ESTIMATE 


In the investigation of J^(12), the quantity 


x»oo >*oo 

$T = J a>j(r) / o> 2 (r J )dr’ dr 
* / R Jx 


was neglected on the grounds that it could be shown to be of order l/R unless both 
2 2 2 

kg — kj and k — kg. Essentially the same integral was neglected in the detailed study 
of Jq for the same reason. In this appendix, a rigorous estimate of the magnitude 
of & will be given. 

If R is large enough, asymptotic forms may be used for and w 2 , and $ji can 
be written as the sum of eight integrals all of the general form 


P = 


[9[ 


j-M s 

MP 


ds dr 


(Al) 


The various values of A and p will be ±ik 2 ± ik and ±ikj - ik, respectively, and thus 
they satisfy the conditions 


#e(p) > 0 
[ft e(p + A) > 0 


(A2) 


It is convenient to write P in the form 


where 


P = 


L 


g(r) 

l+a 


dr 


g(r) = e 


-Ar 


/' 


,-(M+ A)r 


'T 

-(p+ A)r r o -| 

g(r) = [$(1, 1 - 0; pr) - r(l - /3)(pr)' e Mr J 

fir? 


,-Ps 

: ds = 

,1+0 J 


*(1,1 - /3; pr) 


(A3) 


(A4) 


(A5) 
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Equation (A4) may be used to show that 


h(r) = prg(r) (A6) 

is bounded as r — °o if p. is fixed; and equation (A5) may be used to show that g(r) is 
bounded for all r > R even if p — 0. Thus, if p is fixed. 



(Note that this result is independent of the value of A. ) If p — 0, however, it is neces- 
sary to integrate equation (Al) once by parts and write 


,-(X+jLi)r 


+ (1 + a)g(r) 


Now, since g(r) is bounded even for p — 0, the integral is of order l/R and 

(A 9) 

Equations (A7) and (A9) together show that P is of order l/R unless both A and p 
approach zero. 
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